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Spin relaxation from a triplet excited state to a singlet ground state in a semiconductor quan- 
tum dot is studied by employing an electrical pump-and-probe method. Spin relaxation occurs via 
cotunneling when the tunneling rate is relatively large, confirmed by a characteristic square depen- 
dence of the relaxation rate on the tunneling rate. When cotunneling is suppressed by reducing the 
tunneling rate, the intrinsic spin relaxation is dominated by spin-orbit interaction. We discuss a 
selection rule of the spin-orbit interaction based on the observed double-exponential decay of the 
triplet state. 



Electron spin in semiconductors has been a focus of 
research in the context of spintronics, in which spin is 
manipulated with spin-orbit coupling |2| , and of quan- 
tum computation, in which spin carries a quantum in- 
formation |3|. In contrast to two-dimensional electron 
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FIG. 1: (a) SEM picture of the device together with a 
schematic of the measurement set-up. (b) DC current, /, as a 
function of Vl and B measured at V s d =0.15 mV in the first 
cool-down. / for the even N data is divided by 10. Evolution 
of the current peak is sketched around the transition fields 
(triangles), (c) dl/dVh at V s d = 1.2 mV taken at the dashed 
square region shown in (b). Ast scale shows the region where 
the spin relaxation measurement is conducted, (d) dl/dVh at 
V s d = 1.0 mV measured in the second cool-down. The vertical 
dashed lines denote spin transitions of the ground state. A 
blue circle marks an anti-crossing between two singlet states. 
Intensity for the high B data is enhanced by multiplying some 
smooth numerical function. 



gas (2DEG) with continuum density of states, electron 
spin in a quantum dot (QD) is basically free from elas- 
tic scattering, and the resulting long-lived spin states are 
favorable for spin-based applications. Indeed, relaxation 
times of more than 100 [is have been reported in QDs 
between Zeeman sublevels 0, 0, 0, as well as between 
a spin triplet and a singlet state 0, E3|- These relax- 
ation processes have been discussed in terms of either 
spin-orbit interaction or the cotunneling effect. In this 
Letter, we study spin relaxation from a triplet state to 
a singlet state in a lateral QD, in which all the relevant 
parameters can be controlled with the gate voltages. We 
observe smooth transition of the relaxation mechanism 
from the cotunneling regime to the spin-orbit regime by 
changing tunneling rates. The decay of the excited triplet 
state follows a single exponential curve in most condi- 
tions, but double exponential behavior is observed at a 
particular magnetic field where the triplet state crosses 
another state. This might be related to the long-lived 
spin-entanglement state under strong spin-orbit interac- 
tion. 

Figure ^a) shows a scanning electron micrograph 
(SEM) image of our QD device. The AlGaAs/GaAs 
2DEG is constricted by combined dry-etching and sur- 
face Shottky gates. We use only the three gates on the 
right-hand side to form a single QD as shown by the 
white circle. All the measurements are performed in a 
dilution refrigerator at ~90 mK with magnetic field, £>, 
applied perpendicularly to the 2DEG. 

The dot used in this study has charging energy of 
U ~ 2 meV and electron number TV ~ 8. When the 
magnetic field is not very small (B > 0.4 T), electron 
orbitals in a QD can be classified by Landau level (LL) 
index, to which they approach in the high field limit 0|. 
Many-body correction of direct and exchange Coulomb 
interactions induces spin and orbital transitions associ- 
ated with different LLs in low magnetic field, but with the 
same LL in high magnetic field [l£|. Figure 1(b) shows 
an observed DC current, /, through the dot as a func- 
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tion of the left gate voltage, Vl> and B. The two stripes 
show a pair-wise motion with B reflecting spin degener- 
acy. The lower stripe, corresponding to odd TV, involves 
a level crossing (denoted by a solid triangle), associated 
with two orbit als in different LLs. The upper stripe for 
even N involves two spin transitions (open triangles) un- 
der Coulomb interactions. The ground state for even N 
is assigned to be spin triplet between the two transition 
fields, otherwise spin singlet state These spin states 
can be observed in the excitation spectrum of Fig. [He), 
in which the derivative of the current, dl/dV^, with a 
large V s d = 1.2 mV, is plotted as a function of Vl and B. 
Some excited states as well as the ground state that fall 
within the source-drain transport window are observed. 
We study spin relaxation from the triplet excited state 
to the singlet ground state (denoted by the arrow) sep- 
arated by energy, A$t- This relaxation involves orbital 
change between different LLs (inter LL transition). Fig- 
ure 1(d) is another excitation spectrum taken with the 
same device but in the second cool-down. In addition to 
the similar singlet-triplet transitions at B ~ 0.5 T, four 
spin-flip transitions are resolved at B = 2 — 3 T until 
the system enters a stable totally spin-polarized regime 
(y = 1), from which N = 8 is estimated pj]. We also 
study spin triplet-singlet relaxation (denoted by the ar- 
row) that involves orbital change within the same LL 
(intra LL transition). The measurement in the 1st (2nd) 
cool-down with relatively fast (slow) cooling speed re- 
sulted in moderate tunnel rates and spin transition fields 
suitable for studying spin relaxation involving inter (in- 
tra) LL transition, but both sets of data show similar 
characteristics. 

An electrical pump-and-probe measurement is per- 
formed by applying two-step square pulses to the plunger 
gate 7]. First, the singlet and triplet states are emptied 
by lifting both states above the Fermi energy as shown 
in the inset to Fig. Efa) (initialization). The duration 
of this initialization is t\. Next, both states are pulled 
down below the Fermi energy as shown in the inset to 
Fig. Efb). Then, only one electron can enter the dot be- 
cause of the Coulomb blockade. This electron, if it popu- 
lates the triplet state with a probability, P, is allowed to 
relax to the singlet ground state while the pulse height 
is kept at this condition during the wait time, t^. The 
three triplet sublevels with Sz = ±1, are presumably 
populated with an equal probability since the Zeeman 
splitting is negligibly small and the tunneling probabil- 
ity does not depend on Sz- Finally, the pulse height is 
adjusted so that only the triplet state is within the trans- 
port window of 150 jaeV defined by the Fermi energy of 
the left and right leads (read-out). Then, the electron can 
contribute to the current only if it remains in the triplet 
state after t^. This read-out pulse width, t m , is fixed to 
500 ns. Actually, several electrons (1/(1-P)) flow during 
this time for the unrelaxed case. Therefore, the average 
number of tunneling electrons per one pulse cycle, n t , fol- 
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FIG. 2: (a) nt as a function of the initialization time, t\, with 
A ST = 300 ueV (B = 0.6 T in the first cool-down). The inset 
shows the energy diagram of the initialization process, (b) n t 
as a function of the wait time, th- The inset shows the energy 
diagram of the relaxation process. 

lows an exponential decay nt = j^p exp(— th/r s ), from 
which the spin relaxation time, r s , can be determined 0. 

Figure 2(b) shows observed n t as a function of th at 
Ast = 300 ueV. n t shows a single exponential decay with 
r s = 90 us. The relatively large n t (th = 0) ~ 8, corre- 
sponding to P c± 0.89, comes from the fact that an injec- 
tion into the triplet (singlet) state is more (less) effective 
because an electron is added to the outer (inner) orbital 
with a larger (smaller) tunneling rate in this magnetic 
field region |9|. 

We can also determine the total tunneling rate, r^t 
(= Tl + Tr) by changing the initialization pulse width, 
t\. Here, Tl (Tr) is the tunneling rate for the left 
(right) barrier, which is changed by Vl (Vr). Fig- 
ure Efa) shows an example of such measurement. The 
rise time corresponds to the escape time from the sin- 
glet ground state through both tunneling barriers, and 
r to t is estimated by fitting the data to the expression 
n t = j^p{l — exp(— tiTtot)}- r tot is successfully changed 
between 1 x 10 8 and 3x 10 9 s _1 by changing the gate volt- 
ages. 

When r to t is large, higher-order tunneling, or cotun- 
neling, is quite effective in causing an exchange of elec- 
trons having opposite spins between the dot and the lead 
electrodes, resulting in a spin relaxation. According to 
the second-order perturbation theory, the cotunneling 
rate, r c ~J~, is approximately given by 

r~l = AsT^rtJ 2 ^! 1 + 5^) 2 /h, (1) 

where 5- and £ + are energies required to excite the TV 
electron triplet state to TV — 1 and N + 1 electron virtual 
states, respectively pj. rj* ot is the effective tunneling 
rate for the cotunneling process from the triplet to the 
singlet state through either virtual state, while the ex- 
perimentally obtained r to t measures the tunneling rate 
from the singlet state to N — 1 electron state. Figure 
[2a) shows observed Tg as a function of r to t for repre- 
sentative conditions. The data points at large r to t are 
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FIG. 3: (a) Log- log plot of r s as a function of r to t. t s for inter 
LL is measured at B = 0.6 T (Ast = 300 /xeV) in the 1st 
cool-down, while that for intra LL is measured at B — 1.2 T 
(Ast = 380 /xeV) in the 2nd cool-down. The solid lines are 
fitted to the data, (b) The coefficient a for the cotunneling 
component as a function of Ast- The curve is calculated with 
an effective tunneling rate, r* Q t — P^tot- 



almost parallel with the dotted line having a slope —2, 
i.e. r~ ] — cd^ ot , which is consistent with Eq. ([T]). as- 
suming a linear relation T^ ot = /?r to t- Figure 0Jb) com- 
pares the observed Ast dependence of a with calculated 
a = {hAsTiSZ 1 + ^ 1 ) 2 /3 2 }/4tt 2 . We take 5- = 5+ = 
(U/2 — Ast) in the calculation, which approximates the 
experimental conditions. (3 of 0.3 gives a reasonable fit 
to the experimental results. Therefore, spin relaxation in 
the large r tot regime can be well explained by the stan- 
dard cotunneling theory. 

It is seen in Fig. CJa) that, when r to t is reduced, r s 
increases and eventually saturates. In this regime, the 
cotunneling process is suppressed, and inelastic spin re- 
laxation is dominated by phonon emission under the spin- 
orbit coupling effect rather than by coupling to nuclear 
spins, etc 0, 0, The solid lines in Fig. |3fa) are 

the curves (l/r so +l/r cot ) _1 fitted to the data. Here, r so 
is the relaxation time due to the spin-orbit interaction, 
which is independent of T to t • 

Figure Hfb) shows Ast dependence of r s measured in 
the small T to t regime where spin-orbit interaction is dom- 
inant. r s is almost constant in a wide Ast regime (except 
at a dip around Ast ~ 380/ieV) and tends to increase 
when Ast < 200 /ieV for both inter and intra LL data. 
This feature might have arisen from the phonon emission 
spectra in a QD in the presence of spin-orbit coupling. 
The phonon emission rate is maximized when the phonon 
wavelength is comparable to the dot size (phonon energy 
of 300 fieV for the dot size of 30 nm in the crystal growth 
direction) 0,E3- 

The similarity of the Ast dependence 
of the inter and intra LL data supports this crude model. 
Longer r s is observed for the intra LL case than for the 
inter LL case in the whole Ast region explored. This 
might reflect the different orbital quantum numbers in- 
volved in each case, which are relevant to the orbital 
effect on the phonon emission and spin-orbit interaction 



[lilQ- However, the precise mechanism is not clear yet. 

Generally, in a one-electron system, the spin- up (down) 
state of one orbital is coupled with the spin-down (up) 
state of the other orbital when spin-orbit coupling is con- 
sidered between the two orbitals This coupling gives 
rise to finite phonon emission probability between Zee- 
man sublevels. In the case of two-electron singlet /triplet 
states, singlet state (\S)) is coupled with two of the 
triplet sublevels (|T+) and |T_)) having S z = ±1 but 
not with the other sublevel (|T )) having S z = [HE! 
Therefore, relaxation from the triplet to the singlet state 
should have a selection rule in which |Xb) state is still 
free from the spin-orbit relaxation mechanism as shown 
in Fig. E{d). This simple argument applies when the 
singlet-triplet energies are so close to each other that 
coupling with other states is negligible. Unfortunately, 
our pump-and-probe technique is not available for the 
small Ast regime because the minimum energy resolu- 
tion is about 100 fieV. When another singlet excited state 
is involved, however, clear selectivity may appear in the 
vicinity of the level crossing between the triplet and the 
singlet excited states. 

Actually, as shown in Fig.EJc), we do see a level cross- 
ing with an unknown state (X) at about B c± 0.52 T 
(Ast — 400 fieV) that could be a singlet state. When 
the dot potential has no rotational symmetry, large anti- 
crossing is expected between the same spin states [lr^ . 
Typical anti-crossing energy between the same spin states 
in our non-circular dot device is about 150 jaeV (for in- 
stance, marked by the circle in Fig.^d)). The unresolved 
anti-crossing between T and X states (< 100/ieV) im- 
plies that X is a singlet state. As shown by the arrow in 
Fig-Efb), a sharp dip in r s is observed at Ast — 380/ieV 
close to the X-T crossing point. The width of the dip 
is as narrow as about 20 /ieV in Ast (~ 10 mT in 
B) The dip could be attributed to strong spin- 

orbit coupling around the crossing point, resulting in 
the short relaxation time, as theoretically predicted in 
Ref. 14]. It should be noted that the decay of the 
pulse-induced current shows a non-single exponential be- 
havior around the dip as shown in Fig. Bid), while sin- 
gle exponential decay is always observed at other condi- 
tions, e.g., at Ast = 300/ieV shown in Fig. E^c). The 
decay characteristic in Fig. ^d) can be very well fit- 
ted with a double exponential funcion (the solid line); 
C7iexp(-t h /r so ) + C 2 exp(-th/r C ot). Here, Ci, C 2 and 
r so are fitting parameters, and r cot is an input parameter 
(= 810 fis) determined by extrapolating the r to t depen- 
dence of r s in the cotunneling regime to the present value 
of r to t • The fast component of the double exponential de- 
cay can be assigned to the relaxation from |T+) and |T_) 
via spin-orbit coupling, while the slow component to the 
relaxation from \Tq) via higher order spin-orbit coupling 
or the remaining cotunneling contribution. We find that 
the obtained ratio of C\jCi is 2.0 and r so is 60 /is. The 
ratio C\ / C2 obtained at slightly different r to t ranges be- 
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FIG. 4: (a) Ast dependence of the excitation probability, P, 
for the inter LL relaxation (1st cool-down), (b) Ast depen- 
dence of the r s for the inter and intra LL relaxations. Data 
with open symbols are obtained by fitting a single exponen- 
tial function, while the solid circle represents the fast compo- 
nent of the double exponential function. The dotted lines are 
guides for the eye. (c) A logarithmic plot of n t vs. th with 
r to t = 2.0 x 10 8 s" 1 at A ST = 300 /xeV (B = 0.6 T). The solid 
line is an exponential function fitted to the data, (d) A loga- 
rithmic plot of nt vs. th with Ttot = 2.7 x 10 8 s _1 at Ast = 
380 fieV (B = 0.55 T). The solid line is a double-exponential 
function fitted to the data. The dotted and dashed lines 
are the fast and slow components, respectively. The inset 
schematically shows allowed and forbidden transitions from 
the triplet sublevels to the singlet state. 

tween 1.6 and 2.0. These values of C1/C2 are close to 
2, which is expected for an equal population of the three 
triplet sublevels. 

The above observations agree well with the selection 
rule for spin-orbit coupling that is enhanced in the vicin- 
ity of the X-T crossing. However, we cannot safely rule 
out other possibilities like populating an X state in addi- 
tion to the triplet state of interest. Indeed, a very small 
increase (~ 1%) in the excitation probability P is noted 
at around B = 0.55 T as shown by the arrow in Fig.0|a), 
which could be due to injection into the X state. How- 
ever, this effect is too small to explain the observed dou- 
ble exponential behavior. Another possibility might be 
coupling to nuclear spins, which often appears at the level 
coincidence of different spin states . 

In summary, we have studied spin relaxation dynamics 
from the triplet excited state to the singlet ground state 
in a lateral QD. The dominant spin relaxation mech- 
anism is cotunneling at a large tunneling rate, and it 
changes to spin-orbit interaction when cotunneling is sup- 
pressed. The observed double exponential decay char- 
acteristic could reflect the selection rule for the singlet- 
triplet transition mediated by spin-orbit interaction. Fur- 
ther investigation is required to prove this is the case. 



We can think of an "entanglement generator" using 
this selection rule: The singlet ground state |5) = 
IT) a U)a holds a spin pair in an orbital a, while the triplet 
state contains non-entangled states |T+) = |t) a |T)& 
and |T_) = ||) a and entangled state |T ) = 

7f (IT) a \i)b + I Da ID&)> with an electron in each of the 
orbital a and b. At a proper waiting time after the elec- 
tron injection [e.g. th = 300 /is in the case of Fig. 0(d)], 
the system is left in the entangled triplet state |Tq) with 
a probability of 18 % [n t (t h = 300/is) x (1 - P)] or oth- 
erwise in the singlet ground state \S). Our pulse mea- 
surement is based on the extraction of an electron from 
the unrelaxed \Tq) state (the outer orbital with the high 
probability P), and thus this scheme can be used to gen- 
erate or analyze an entangled spin pair by detecting the 
extracted electron with a sensitive electrometer |g, |S| • 
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